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^ ' Abstract 

i^ , Pedestal profiles that span the ELM cycle have been obtained and used to test 

the idea that the pedestal pressure gradient in MAST is limited by the onset of Ki- 
netic Ballooning Modes (KBMs). During the inter-ELM period of a regularly type 
^ . I ELM-ing discharge on MAST, the pressure pedestal height and width increase 

^\ \ together while the pressure gradient increases by only 15 % during the ELM cycle. 

Stability analyses show that the pedestal region over which infinite-n ballooning 



^r^ . modes are unstable also broadens during the ELM cycle. To test the relationship 



between the width of the region that is unstable to n = cxd ideal magnetohydro- 
C^ ' dynamic ballooning modes and KBMs the gyrokinetic code, GS2, has been used 

for microstability analysis of the edge plasma region in MAST. The gyrokinetic 
simulations find that KBM modes with twisting parity are the dominant microin- 
stabilities in the steep pedestal region, with a transition to tearing parity modes 
^ ■ in the shallower pressure gradient core region immediately inside the pedestal top. 

The region over which KBMs are unstable increases during the ELM cycle, and a 
good correlation is found between the region where KBMs dominate and the region 
that is unstable to infinite-n ideal ballooning modes. 



1 Introduction 

The reference scenario for ITER [1] is the high confinement mode (or H-mode) [2]. The 
improved confinement is the result of a narrow transport barrier that forms at the plasma 
edge. Due to the strong coupling between pedestal and core confinement, which has been 
observed in many devices [3], the fusion performance of ITER will strongly depend on the 
achievable pedestal pressure. Therefore, an accurate prediction of the pedestal height is 
essential for the prediction and optimization of ITER performance. A class of instabilities 
called Edge Localized Modes or ELMs [1], [5] limits the maximum pressure pedestal 
height, Ppcd; that can be achieved for a given pedestal width. Hence, understanding the 
width of the barrier is very important in determining the maximum pressure that can 
be achieved. An MHD stability analysis has been used to show that in order for ITER 
to attain a temperature pedestal of 4keV, which is thought to be required for ITER to 
attain its fusion power output objective, the pedestal width would have to be 2.5% of 
the minor radius [6], [7]. Several empirical scalings for the pedestal width (Aped) have 
been proposed based on dimensionless parameters such as the normalized poloidal ion 
Larmor radius (PpoJ or pressure (/3poi), but a strong co-linearity between pedestal values 
of /3poi and p^^y limits the power of these datasets to discriminate between very different 
scalings: e.g. early data from DIII-D could be fit equally well with either Aped = (Ppoi)^'^^ 
or Aped = (/SpoO^'"^ [S]- Experiments on JT-60U using hydrogen and deuterium plasmas 
found a scahng of the form Aped oc (Ppoi)'''^(/5poi)°'^ [9]. Recent similarity experiments on 
DIII-D, ASDEX Upgrade and JET also obtain a pedestal width scaling, in normalised 
poloidal flux, of the form A*^ oc f3'^'^i [10]. Previous studies on MAST have shown no 
evidence for an increase of the width with p* [ but a clear increase of width with /3poi [H] . 
Fits of the distributions of MAST pedestal width measurements to the form Aped = ^/^^i 
yielded B ^ 0.5. A similar scaling of the pedestal width has also been observed on DIII-D, 
C-MOD and AUG [12]. The scaling Aped oc (3^^^ has been used to provide basic input to 
the EPED model for pedestal evolution [13]. The EPED model proposes that: drift wave 
turbulence is suppressed in the pedestal region by sheared flow; kinetic ballooning mode 
(KBM) turbulence constrains the pedestal to a critical normalized pressure gradient; and 
an ELM is triggered by the onset of the intermediate n (~ 10 — 20) peeling-ballooning 
mode instability. This model provides a prediction for the pressure pedestal width and 
height just before the ELM crash. EPED was successfully compared to pedestal conditions 
in DIII-D and a range of other tokamaks |T2],[13]. The model has since been extended 
to exploit the fact that KBMs are predicted to have an onset condition similar to that of 
ideal infinite-n ballooning modes [H] , where n is the toroidal mode number. In the more 
recent version of EPED, the infinite-n ideal MHD ballooning mode code, BALOO [T5] . 
has been used to compute the threshold pressure gradient set by the condition that either 
1% of the pressure profile in normalized poloidal flux [H] or 50% of the steep pedestal 
gradient region [TB] is unstable to infinite-n ideal ballooning modes. This technique has 
found good agreement between the predicted and observed pedestal widths and heights 
obtained prior to type I ELMs on DIII-D [H] and JET [17]. 

In this paper, the ideas incorporated in the EPED model have been tested using data 
from the spherical tokamak MAST. In order to construct a pedestal model for MAST, 
the evolution of the pedestal profiles has been studied as a function of time during the 



inter-ELM period, and this is described in Section 2. In Section 3 the HELENA code [H 
has been used to reconstruct equihbria during the ELM cycle and to determine infinite-n 
ideal ballooning stability, and the ELITE code [T9l |20] has been used to calculate the 
finite-n stability boundary. In Section 4 the local gyrokinetic code, GS2 [21], has been 
used to perform linear microstability analysis in the pedestal region and to test whether 
the region unstable to infinite-n ballooning modes is related to the region where KBMs 
are unstable in MAST (i.e. to test whether ideal infinite-n MHD ballooning analysis can 
be used as a proxy for KBM stability in MAST). 

2 Inter-ELM pressure gradient evolution on MAST 

The data examined in this paper were obtained from a set of three MAST discharges 
^24459, ^24452 and #24763. All discharges were in a double null configuration with 
a plasma current of 850 kA, a toroidal magnetic field on axis of 0.585 T and heated by 
3.4 MW of neutral beam injection. The timeslices examined were obtained during the 
inter-ELM period of a regularly type I ELM-ing discharge. 

The data used in this paper were taken from a Thomson scattering system with 130 
spatial points [22l [23] . This system combines eight Nd:YAG lasers and measures Te and 
Ug profiles approximately every 4.2ms. This sampling rate is not normally sufficient to 
diagnose an individual ELM period (6-12ms), but arranging the profiles obtained over a 
number of plasmas as a function of time with respect to the ELM cycle allowed the general 
cycle of periodic ELMs to be diagnosed. Each spatial point in the Thomson scattering 
system has a scattering length of 10mm. At the Low Field Side (LPS) edge of the plasma 
this results in 10mm radial resolution. An, as the laser path is in an approximately radial 
direction. Close to the High Field Side (HFS) edge of the plasma, the laser is at an angle 
to the radial direction and this results in much higher (3-5mm) radial resolution. 

Measured edge profiles were fitted to an analytic parameterisation, P{r; a), convolved 
with the instrument function to account for the spatial resolution of the diagnostic. P{r; a) 
is based on a modified hyperbolic tangent function, mtanh, and is defined 



P(r; a) 
where 
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using as fit parameters: pedestal height, Oped; scrape off layer height, Osoi; transport 
barrier position, adh] width, ca; and core slope, Csiope- The pedestal width for profile 
quantity a, A^, is determined from these parameters as Aq = 4a a- Electron temperature 
and density pedestal width measurements on the LFS were typically close to the limits of 
spatial resolution, giving Axe.ne/Aij ~ 1 — 2: but on the HFS there were sufficient data 
points in the pedestal region to reduce the uncertainties, giving ATe,ne/^R ~ 3 — 4. All 
data examined in this paper were therefore obtained from measurements at the HFS of 
MAST. 



Three profiles obtained during a single ELM cycle are shown in figure [H The Te(r) 
and ne(r) profiles were first fitted in real space using the parameterisation of equation ([1]) 
convolved with the instrument response function. These fitted profiles were then mapped 
to normalised poloidal fiux space {"^n), where a further fit was used to obtain pressure 
profile fit parameters in \E'Af- The electron pressure gradient, dPe/d'^N, was computed 
from this fit. 

During the ELM cycle examined in figure [H the change in the temperature profile is 
small. The density profile shows a large change, with Helped and A„e increasing through- 
out the ELM cycle. The resulting pressure pedestal change, dominated by the growing 
density pedestal, shows increasing Pe,ped and expanding Ape. The pressure gradient pro- 
file evolution during the ELM cycle shows an expanding region of high pressure gradient 
(with little change in the peak pressure gradient), indicating that a transport barrier is 
expanding inwards towards the core. 

Examining the peak pressure gradient obtained in this manner over several inter-ELM 
periods, and plotting the result as a function of time during the ELM cycle (exploiting the 
fact that the ELMs are periodic), results in the data shown in figure |2j From this graph 
we see that immediately after the ELM, the drop in pressure gradient recovers rapidly 
during the first 10% of the ELM cycle. During the remainder of the cycle there is only a 
weak upwards trend in the peak pressure gradient. 

A more complete picture of the pedestal evolution in the inter-ELM period was ob- 
tained by combining the 50 profiles from the three MAST discharges. For each timeslice, 
the mtanh fit parameters for electron density, electron temperature and electron pressure 
profiles in \E'7v space, have been obtained as previously described. Least squares linear 
fits were found to represent well the time evolution of each mtanh parameter over the 
ELM cycle. Such linear fits were used to evaluate the mtanh profile parameters at five 
normalised times during the ELM cycle, t = (0.1,0.3,0.5,0.7,0.9). Table [T] gives the in- 
terpolated values of the mtanh fit parameters for the electron pressure profile at these 
times, and figure [3] shows the reconstructed electron profiles. The peak electron pressure 
gradient from each of these profiles is compared with the peak electron pressure gradients 
determined directly from mtanh fits to individual profiles in figure [2J 

Figures El^a) and (b) show little change in the electron temperature pedestal height, 
whilst Atc, A„e and the density pedestal height increase during the ELM cycle. It can be 
seen from figure [3]^c) that the electron pressure pedestal height and width increase during 
the ELM cycle. Equilibrium reconstruction shows little change in the fiux expansion at 
the HFS during the ELM cycle, and therefore the Ape expansion in real space tracks 
that in "^^ space. The calculated electron pressure gradient is shown in figure IHl^d). The 
increase in Pe,ped over the ELM cycle is mainly associated with the expansion in Ape, and 
only partly with a modest increase in the pressure gradient, which increases by < 20% 
over the ELM cycle. The high pressure gradient pedestal region expands during the ELM 
cycle: the innermost surface where dp/d'^N > lOkPa moves inwards from "^j^ = 0.975 at 
t = 0.1 to \E'7v = 0.95 at t = 0.9, while the location of the peak pressure gradient also 
moves inwards from "^j^ = 0.988 at t = 0.1, to ^jy = 0.978 at t = 0.9 (where t denotes 
normalised time during the ELM cycle). Similar observations were made during regular 
type I ELM-ing periods in other plasma scenarios on MAST, DIII-D [13] and AUG p5]. 



3 MHD Stability Limits During the ELM Cycle 

Measured electron temperature and density profiles were used to reconstruct MHD equi- 
libria to describe the edge plasma with high precision. Five such equilibria were produced 
to span the ELM cycle. Pressure profiles were obtained from the rie and Tg profiles of 
figure 3 and model profiles for ion temperature and density, as ion measurements were 
not available with sufficient temporal and spatial resolution. A uniform effective charge 
state, Zeff = 2, [^ was assumed, for a plasma with Carbon as the main impurity species so 
that rii = rieiy — Zeff)/Q. Two ion temperature models were adopted: Tj = Tg across the 
whole plasma including the pedestal region; and setting Tj = Tg in the core, but extrapo- 
lating Ti from the pedestal top to the separatrix using the electron temperature gradient 
from the core side of the pedestal top, to avoid a steep edge gradient in Tj. The plasma 
boundary shape, total plasma current and vacuum magnetic field were taken from MAST 
equilibrium reconstructions using EFIT |26]. The inductively driven current profile was 
assumed to be fully relaxed, and the bootstrap current, which dominates at the edge, was 
calculated using the formula by Sauter et al. [27], [28]. These assumptions have allowed 
the reconstruction of high precision equilibria using the HELENA code [18]. These equi- 
libria exploit the available constraints from the TS measurements self-consistently, and 
were used for the MHD and gyrokinetic stabihty analyses presented here and in Section HI 
respectively. 

In order to determine the limiting pressure gradients for ideal MHD instability, scans 
in the equilibrium pedestal pressure gradient have been performed, either by scaling the 
temperature pedestal height and gradient (at constant pedestal width) keeping the density 
profile fixed, or by scaling the density pedestal height and gradient (at constant pedestal 
width) keeping the temperature profile fixed. Both approaches gave similar results, but 
here results from the latter approach are presented, as this is closer to the experimental 
situation. The equilibrium scan parameter is the maximum value, amax? of the normalised 
pressure gradient, a, given by [29] : 

-2dV/d^ { V V''^ dp 

where p is the pressure, V is the volume enclosed by the flux surface, Rq is the geometric 
centre of the plasma and ip is the poloidal flux. In scans the pedestal pressure gradient 
was varied by up to a factor 2.5 (while the equilibrium stored energy varied by < 30%), 
and the bootstrap current was always calculated self-consistently for each equilibrium 
using the appropriately modified density and temperature profiles. 

Ideal MHD stability analysis was performed using two codes, depending on the toroidal 
mode number n. The HELENA code used for equilibrium reconstruction, was also used 
to determine infinite-n ideal ballooning stability on each flux surface [H]. The ELITE 
code [19], [20] was used to analyse stability of intermediate n modes. In this study the 
finite n modes were restricted to a maximum value oin = 25, which roughly corresponds 
to the maximum n observed in experimental measurements of ELM filaments in MAST 
The lowest n considered was ra = 5. 



^Zcff = 2 is consistent with typical experimental values in the MAST pedestal. The core Z^s is 
typically in the range 1-1.5. 



All five experimental equilibria are unstable to infinite-n ideal ballooning modes. The 
radial region that is unstable to n = oo modes is plotted in figure H] as a function of 
normalised time in the ELM cycle, and will be compared to the region unstable to KBMs 
in Section HI Infinite-n ballooning stability was not found to be sensitive to the choice of 
ion temperature model in the pedestal region. All other MHD and gyrokinetic stability 
analyses were also insensitive to the ion temperature model, so the remaining stability 
results presented in this paper assumed Tj = Tg throughout the plasma. 

Figure [5] shows the growth rate of the most unstable intermediate n ideal MHD mode, 
n = 25, as a function of Omax, for scans around the first and last equilibria in the ELM 
cycle. The threshold for the finite-n instability, Oc, corresponds to the Omax value above 
which the mode becomes unstable with a rapidly rising growth rate. Figure [5] illustrates 
that while the experimental pedestal pressure gradient, aexp, is far below ac at the start 
of the ELM cycle, it is slightly above «c by the end of the cycle. 

Figure 3 shows that the maximum pressure gradient increases only slightly during the 
ELM cycle: the maximum of the edge bootstrap current density (oc Vp) also changes 
very little. The main change during the ELM cycle is the broadening of the pedestal. 
Figure |6] plots Ooxp and ac as functions of normalised time during the ELM cycle. While 
the experimental pressure gradient changes very little, the stability limit is lowered sig- 
nificantly as the pedestal broadens. This can be understood in terms of stabilising finite 
n corrections in ballooning theory [31], which are larger for the equilibria with narrower 
pedestals from earlier in the ELM cycle. It has frequently been found useful to present the 
proximity of equilibria to edge localised MHD stability boundaries using plots of MHD 
stability versus edge pressure gradient and current density (a^ax, j^i) [22] • Such diagrams 
have considerable value, but their description of stability dynamics over the full ELM 
cycle is limited by the difficulty of capturing the full equilibrium evolution in a simple 
2-parameter model. 

Figure [7] shows the experimental pressure pedestal height as a function of time during 
the ELM cycle. Superimposed on this plot are the stability limits derived from the 
infinite-n ideal ballooning and finite-n stability analyses. The finite-ra stability limit is 
well above the experimental value at the start of the ELM cycle, but later it converges on 
the experimental value just before an ELM, which is consistent with the theory that type 
I ELMs are triggered by an ideal finite-n instability. A stability criterion, based on 100% 
of the pedestal regioncl being unstable to infinite-n ideal MHD ballooning modes, tracks 
the experimental pedestal pressure height to within 20% throughout the MAST ELM 
cycle. Pedestal pressure gradient limits derived from criteria that the region unstable to 
n = oc ballooning modes be wider than (i) 1% of the total poloidal fiux, or (ii) 50% of 
the pedestal, however, lie well below the experimental pedestal pressure throughout the 
cycle. 

In addition to analysing the stability of equilibria during the ELM cycle, an analysis 
similar to that used in the EPED model [13] was repeated. Edge stability was determined 
for a range of equilibria, using a scan around the reference equilibrium just prior to the 
ELM crash. In this scan the density and temperature pedestal widths and the density 
pedestal height were scaled together around their measured values, keeping the density 



^Thc width of the pedestal region is defined: Apod = (A„^ + ArJ/2. 



gradient constant. Then, for each pedestal width Aped, the temperature pedestal height 
was varied to determine stability boundaries. Figure |8] plots various stability limits as 
functions of the scan parameter Aped: the finite-n limit; and limits based on 1% of the 
poloidal flux, and 50% and 100% of the pedestal, respectively, being unstable to infinite-n 
ideal ballooning modes. The pedestal top pressure limit due to finite n modes increases 
with pedestal width, but more slowly than the scanned pedestal top pressure (which is 
oc Aped if the pressure gradient is constant). The limit requiring 100% of the pedestal 
to be unstable to ra = oo ideal ballooning modes, crosses the finite n limit close to the 
experimental pedestal width just before an ELM. While this could in principle provide an 
EPED-like model for the pedestal height and width in MAST |j, the figure of 100% has 
no theoretical justification. It will be necessary to understand this parameter to make 
predictions for future devices. 

4 Microstability analysis during the ELM cycle 

Kinetic descriptions can improve on MHD models of edge instabilities, and the gyroki- 
netic model is exploited here. The H-mode pedestal is associated with a strong pressure 
gradient that can drive electromagnetic instabilities, like the KBM, where the magnetic 
perturbations are crucial to the instability mechanism. Experimental measurements sug- 
gest that the KBM may play a very important role in the H-mode edge, by setting a hard 
limit on the maximum pressure gradient achievable in the pedestal region [H]. Early 
work has shown that collisionless KBMs and ideal MHD infinite-n ballooning modes are 
described by related equations, and that these modes have similar character [33]. This 
suggests that it can be reasonable to use the stability of the infinite-n ballooning mode as 
a proxy for KBM stability. Nevertheless, finite Larmor radius (FLR) effects, only included 
in the kinetic treatment, are often stabilising [3l], and the region unstable to infinite-n 
ideal MHD ballooning modes can be significantly broader than the region unstable to 
KBMs. The importance of FLR effects can be measured by the parameter A: 

K=^Ml^ (3) 

2Lp ^ ' 

where Lp is the pressure gradient scale length, ky is the perturbation wavenumber in the 
flux surface and perpendicular to magnetic field lines, and pi is the ion Larmor radius. 
When A <^ 1 (corresponding to Lp ^ {kyPi)'^R), FLR effects are negligible and the KBM 
stability limit is close to that of the ideal infinite-n ballooning mode: when A ~ 0(1), 
however, FLR effects become significant and should not be neglected [31] • This suggests 
that the ideal infinite-n calculation can only describe KBM stability consistently in regions 
of shallow pressure gradient, and raises questions as to the validity of using infinite-n ideal 
ballooning stability as a proxy for KBM stability in the pedestal region. 

In order to investigate this, a local electromagnetic gyrokinetic analysis has been used 
to obtain the microstability properties of the MAST edge plasma during the ELM cycle. 



■^A pedestal evolution model cannot be constructed for MAST using a criterion based on either 1% of 
poloidal flux [T3] or 50% of the pedestal [IF] being unstable to infinitc-n ideal ballooning modes, as the 
experimental pressure profile appears to be well above this limit. 



The gyrokinetic model is based on a first order expansion in the small parameter 5 = pi/ L 
(where L is a typical equilibrium gradient scalelength) , which includes FLR corrections 
that are missing from ideal MHD. Clearly gyrokinetics should be less accurate in the 
pedestal than in the core, as the pedestal value of 5 is larger owing to shorter equilibrium 
gradient scale lengths. Using the MAST data at the time corresponding to the fully 
developed pedestal (i.e. just prior to the ELM crash), we obtain 5 ~ 0.3 (ie 5 < 1) at 
'^ N = 0.98, which corresponds to the surface at the mid-radius in the fully developed 
pedestal. At this location it is interesting to note that A ~ 1 for kypi ~ 0.2, so FLR 
corrections should be important for kyPi > 0.2. While conditions in the pedestal limit 
the accuracy of the gyrokinetic model, gyrokinetics is more appropriate than MHD as it 
includes FLR corrections. 

The local gyrokinetic code, GS2, has been used to perform linear microstability anal- 
ysis around the pedestal during the MAST ELM cycle, based on the data described in 
Section 12] and the equilibria prepared for MHD stability analysis in Section 3. GS2 is an 
initial value code, and finds, for specified values of kyPi, the fastest growing microinsta- 
bility and provides its growth rate, 7. GS2 was used to obtain the dominant growth rate 
spectra, which were captured in the range 0.07 < kyPi < 5.5, at 12 radial locations between 
^N = 0.94 and ^n = 0.995, for 5 times during the ELM cycle, t = (0.1, 0.3, 0.5, 0.7, 0.9). 
The characteristics of the dominant mode reveal the radial region where KBMs are ex- 
pected to dominate, and this can be compared directly with the region found to be 
unstable to infinite-n ideal ballooning modes using HELENA. 

Figure M^a) shows the parity factor, Cpar, of the dominant mode in the ky range 
0.07 < kyPi < 5.5, and figure [9]^b) shows the growth rate of the fastest growing mode at 
kyPi = 0.218, which corresponds to the wavenumber of the dominant KBM instability over 
the ELM cycle; both quantities are illustrated as functions of radius and time through 
the ELM cycle. The parity factor is defined: 

\fA\\de\ 

r _ 1 _ 1-^ II I 
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where A\\ is the parallel magnetic vector potential. The dominant modes are found to 
be KBMs with twisting parity (Cpar = 1) in the steep pressure gradient region closer to 
the edge, with a transition to electron temperature gradient driven tearing parity modes 
(Cpar = 0) in the shallower gradient core region inside the top of the pedestal. The switch 
in mode parity is illustrated in figure [TD], which presents the dominant eigenfunctions of 
the electrostatic potential, $, and A\\, on the \I/Ar = 0.98 and '^ jq = 0.95 surfaces at 
t = 0.5. Microtearing modes have previously been reported to be unstable at mid-radius 
in MAST plasmas [35], [36], close to the edge in AUG plasmas [37], and in a conceptual 
fusion power plant based on the spherical tokamak |38j . 

The infinite-n ideal ballooning unstable region (from figure H]) is also indicated on 
figure ini^a), and this corresponds closely to the region where KBMs are the dominant 
modes. This suggests that it is reasonable here to use infinite-n ballooning calculations as 
a proxy for the dominance of KBMs, even though A is finite. The KBM unstable region 
broadens and expands inwards with the strong pressure gradient, as the pedestal width 
increases through the ELM cycle. Figure [Hj^b) shows that the KBM growth rate peaks in 
the steep gradient region close to the top of the pedestal. This peak moves inwards as 

8 



the pedestal expands, and the KBM growth rate reduces at this location after the transit 
of the peak. 

5 Conclusions 

Pedestal profiles of electron temperature and electron density, obtained as a function of 
time during the MAST ELM cycle, have been used to test whether the pressure pedestal 
gradient in MAST is limited by the onset of Kinetic Ballooning Modes (KBMs). In 
this set of discharges, the electron density and electron pressure pedestal heights and 
widths increase during the inter-ELM period, while their gradients remain approximately 
constant. Finite-n stability analyses using the ELITE code show that the finite-n stability 
limited normalised pressure gradient, Uc, lies above the experimental value, Oexp, initially 
just after an ELM, and then falls towards it just before the next ELM: i.e. rather than 
the experimental pressure gradient increasing towards a fixed stability boundary it is the 
pressure gradient stability limit that moves towards the experimental pressure gradient 
due to the increasing pedestal width. The region over which infinite-n modes are unstable 
also broadens during the ELM cycle. Ideal MHD stability boundary calculations track 
the evolution of the pedestal height and width well, provided the limiting criterion used 
is that 100% of the pedestal width be unstable to infinite-n modes. 

The gyrokinetic code, GS2, was used to test whether the MAST pedestal region unsta- 
ble to KBMs corresponds closely to the region that is unstable to infinite-n ideal ballooning 
modes. KBM modes with twisting parity were found to be the dominant microinstabilities 
in the steep pedestal region, with a transition to tearing parity modes in the shallower 
pressure gradient core region immediately inside the pedestal top. The region over which 
KBMs dominate increases during the ELM cycle, and closely corresponds to the region 
unstable to infinite-n ideal ballooning modes. The impact of sheared flow on stability has 
not been considered here, and an assessment of its possible impact is underway. 
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t 



e,etb 



(^ 



'N, 



100 Pe,A {^Pn) Pe,ped (Pa) Pe,slope (Pa) Pc,sol (Pa) 



0.1 0.988±0.0009 0.527 ±0.027 557.5±40.9 0.167±0.010 0.56±0.31 



0.3 0.986±0.0008 

0.5 0.983±0.0009 

0.7 0.981±0.0011 

0.9 0.979±0.0014 



0.606 ±0.026 
0.686 ±0.029 
0.766 ±0.035 
0.846 ± 0.042 



682.5±39.1 
807.6±43.4 
932.7±52.3 
1058.±63.9 



0.166±0.009 
0.164±0.011 
0.162±0.013 
0.160±0.016 



0.47±0.29 
0.38±0.30 
0.29±0.34 
0.20±0.40 



Table 1: The five intanh fit parameters for electron pressure at t = (0.1, 0.3, 0.5, 0.7, 0.9), 
with uncertainties from the least squares linear fits to the evolution of these parameters 
during the inter-ELM period. 
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Figure 1: Profiles obtained during a single ELM cycle, showing the observed electron tem- 
perature and electron density data points in the pedestal, fitted profiles and the inferred 
electron pressure and electron pressure gradient profiles. 
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Figure 2: Evolution of the peak electron pressure gradient during the ELM cycle. The 
circles show the results from f«50 profile measurements. The diamonds show the peak 
pressure gradient from the fitted profiles shown in figure [31 
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Figure 3: Evolution over the ELM cycle of: (a) electron temperature, (b) electron density, 
(c) electron pressure and (d) electron pressure gradient. Profiles were obtained from the 
mtanh fit parameters evaluated at the following normalised times during the ELM cycle: 
t= (0.1,0.3,0.5,0.7,0.9). 
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Figure 4: The unstable region for n = oo ballooning modes as a function of the normalised 
relative time in the ELM cycle for assumptions Tj = Tg (blue squares, solid line) and 



VT, 



i,pcd 



VTg^corc (red circles, dashed line). 
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Figure 5: The growth rate of the most unstable mode, n = 25, as a function of the 
maximum value of the normalised pressure gradient, amax for equilibria from the first 
(solid blue) and last (dashed red) intervals of the ELM cycle. The circle and the star 
show the experimental equilibrium values of amax- Growth rates are normalised to the 
Alfven frequency in the centre of the plasma (ua)- 
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Figure 6: The maximuni normalised pressure gradient in the experimental pedestal, a. 

and the finite n stability limit, a^ 

cycle. 



exp; 



as functions of the normalised time through the ELM 
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Figure 7: The pedestal top pressure of the experiment (black) and for profiles limited 
by the finite n stability limit (blue), and the infinite-n ideal ballooning mode limit using 
three criteria: 1% of the total poloidal flux is unstable (solid red), 50% of the pedestal 
width is unstable (dashed red) and 100% of the pedestal width is unstable (dotted red). 
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Figure 8: The pedestal top pressure as a function of the pedestal width for profiles limited 
by the finite-n stability limit (blue), and the infinite-n ideal ballooning mode limit using 
three criteria: 1% of the poloidal flux is unstable (solid red), 50% of the pedestal width 
is unstable (dashed red), and 100% of the pedestal width is unstable (dotted red). 
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Figure 9: (a) Parity factor, Cpar, of the dominant mode in range 0.07 < kyPi < 5.5, and 
(b) growth rate of the fastest growing mode at kyPi = 0.218, are both shown as functions 
of \l/7v and normahsed time through the ELM cycle. In (a) the infinite-n ideal ballooning 
unstable region is illustrated by black lines and diamonds, and the boundary between the 
regions where twisting and tearing parity modes dominate is indicated by the dashed line. 
In (b) 7 is normahsed to vu/a, where vu is the local ion thermal velocity and a is the 
minor radius. 
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Figure 10: Normalised real part of the $ and A\\ eigenfunctions for the dominant mode 
at the normalised ELM cycle time of 0.5 on surfaces "^n = 0.98 and "^n = 0.95. $ and 
A\\ eigenfunctions are plotted as functions of ballooning angle, 6, and are normalised to 
'^(6'max) and $(6'niax)/'ytj, respectively, where |$| peaks at 6^ = ^max- On the outer surface, 
"^N = 0.98, the dominant mode has kyPi = 0.218, and (a) $ and (b) A\\ exhibit twisting 
parity. On the inner surface, "^n = 0.95, the dominant mode has kyPi = 3.28, and (c) $ 
and (d) A\\ have tearing parity. 
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